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Abstract 

We study general quantum integrable Hamiltonians linear in a coupling constant and represented by finite 
N X N real symmetric matrices. The restriction on the coupling dependence leads to a natural notion of 
nontrivial integrals of motion and classification of integrable families into Types according to the number 
of such integrals. A Type M family in our definition is formed by — AI nontrivial mutually commuting 
operators linear in the coupling. Working from this definition alone, we parameterize Type M operators, 
i.e. resolve the commutation relations, and obtain an exact solution for their eigenvalues and eigenvectors. 
We show that our parameterization covers all Type 1,2, and 3 integrable models and discuss the extent to 
which it is complete for other types. We also present robust numerical observation on the number of energy 
level crossings in Type M integrable systems and analyze the taxonomy of types in the Id Hubbard model. 
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I. INTRODUCTION 



Quantum integrability is usually defined as the existence of a sufficient number of nontriv- 
ial operators commuting with the Hamiltonian. Such operators are interchangeably referred to 
as "conserved currents", "families of commuting Hamiltonians" , integrals of motion, "dynamical 
conservation laws" , or simply "dynamical symmetries" , and are distinct from the ordinary space- 
time or internal space symmetries, e.g. momentum and spin conservation, particle-hole symmetry 
etc. Quantum integrable systems are believed to have a range of properties, which are sometimes 
proposed as alternative notions or tests of integrability - an exact solution for their eigenstates, 
Poissonian energy level statistics, level crossings violating the Wigner-von Neumann non-crossing 
rule in parameter dependent systems, and nondiffractive scattering, see e.g. Ref. Q for a review. 

There are two problems with the above definition. First, it has proved difficult to precisely 
formulate what constitutes a nontrivial commuting operator. For instance, any Hamiltonian com- 
mutes with projectors onto its eigenstates, so the mere statement that there exist commuting 
operators is not meaningful. In classical mechanics one requires functional independence of the 
integrals of motion ^j. This criterion is not similarly useful in quantum mechanics where any two 
commuting matrices are functionally dependent unless there are degeneracies in the spectra of 
both of themy]. Second, in a system with a finite Hilbert space we clearly cannot have an infinite 
number of independent integrals, so we need to specify how many conserved currents are needed 
to say that the system is integrable. This is especially problematic in systems with no well-defined 
classical limit, e.g. the Hubbard model on a finite lattice 4,2], as, unlike the classical mechanics, 
there is no natural notion of the number of degrees of freedom. There is a similar lack of clar- 
ity in delineating the properties of quantum integrable systems. For example, there are poorly 



61] as well as no con- 



understood exceptions to the energy level statistics and level crossing tests 
vincing derivation of the exact solution and other properties from the notion of integrability. This 
is in sharp contrast to the situation in classical mechanics where the Liouville- Arnold theorem Q] 
demonstrates an exact solution by quadratures and quasi-periodic motion on invariant tori to be 
direct consequences of integrability. 

In the present paper we show that most of these problems can be resolved if one considers 
Hamiltonians and conserved currents that depend on a real parameter in a certain fixed way. This 
requirement alone leads to a well-defined notion of quantum integrability allowing us to classify 
and explicitly construct various integrable models, obtain their exact spectra, as well as derive 
other properties, directly from their integrability. Examples of such a real parameter, call it u, in 



condensed matter integrable systems are: Coulomb interaction in the Hubbard model, magnetic 
field in Gaudin magnets with a boundary termQ], anisotropy in the XX Z Heisenberg model, in- 



teraction strength in the BCS modeljg] etc. Dynamical symmetries depend on the same parameter 
in contrast to ordinary symmetries which are n- independent. Cases with no such parameter depen- 
dence, e.g. isotropic Heisenberg magnet and Gaudin model in zero fieldQ], often correspond to a 
more general integrable model at a specific value of the parameter. In the Gaudin and BCS models 
the Hamiltonian and all dynamical symmetries are linear in uj?!. [lol. while in the Hubbard and 
XX Z models the Hamiltonian and the first conserved current are linear, while higher currents are 
higher order polynomials in u 
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Note that in each of the above examples there is a subset of mutually commuting operators, 
including the Hamiltonian, linear in the parameter. Let us therefore consider N x N Hermitian 
operators of the form H{u) = T + uV and require the existence of dynamical symmetries linear in 
u. A simple, but important observation is that this restriction on u-dependence leads to a natural 
notion of a nontrivial commuting partner. A typical ( "non-integrable" ) H{u) will commute only 
with a single (trivial) operator linear in u, namely with H{u) = {a + bu)I -\- cH (u) , where a, 6, and 
c are real numbers and / is an identity operator. The requirement that there exist a nontrivial 
commuting partner linear in u imposes severe restrictions on H[u). Indeed, as we will see below, 
a real symmetric H{u) for which such a partner exists is uniquely specified by a choice of no more 
than AN real parameters, while for a generic H{u) one needs A^(A^ -|- 3)/2 real parameters to fix 
the matrix elements of T and V Thus, Hermitian operators with fixed parameter-dependence 
separate into two distinct classes - those that have nontrivial commuting partners and those that 
do not. 

Further, there is a natural classification of families of parameter-dependent commuting operators 
according to the number K of independent members they contain. Families with the maximum 
possible number K = N (see below) of such operators we term Type 1, with one less than the 
maximum - Type 2, etc. The main results of this paper are as follows. First, we construct Type 
M families of commuting operators W{u) = T* + uV^ for arbitrary M, i.e. we solve commutation 
relations [H^{u), (u)] = for matrix elements of H^{u). Second, we obtain exact eigenvalues 
and eigenstates of each H^[u). We thus derive the exact solution from the integrability. Third, we 
observe that the energy levels of W{u) frequently cross and determine the number of level crossings 
as a function of N and M. Type 1 (maximal) operators were previously constructed in Refs. 1^ and 



19l . Here we analytically obtain all Type 2 and 3 and a sub-class of Type M real symmetric N x N 



operators for arbitrary M > 3. In our construction a Type M commuting family is parameterized 



by 2N + M — 2 real parameters, while, as we show, a generic Type M involves 2N + 2M — 5 
parameters. 

There are typically w-independent symmetries common to the Hamiltonian and its dynamical 
symmetries, e.g. total momentum, particle number, z component of the total spin, etc. It makes 
sense to "factor" all these out, i.e. to go to blocks that correspond to a complete set of quan- 
tum numbers. Matrices we construct essentially represent such "irreducible" blocks of integrable 
Hamiltonians and their dynamical symmetries. As an example, we consider irreducible blocks for 
the Hubbard model with six sites, three spin- up and three spin-down electrons and determine their 
types. 

A toy version of the above approach to integrability based on studying parameter-dependent 
commuting matrices was discussed in Ref. 



2d in the context of level crossings in the Hubbard model. 
A major breakthrough came with Shastry's work [l^ where he constructed Type 1 matrices. The 
advantages of this approach are conceptual simplicity, clear meaning of independent integrals and 
freedom from specific details of particular integrable systems. A significant disadvantage is that 
the H^{u) correspond to isolated sectors/blocks of the Hamiltonian and dynamical symmetries and 
it is difficult to establish a correspondence between matrices H^{u) and underlying operators whose 
blocks they represent. Motivated by Shastry's work, we proposed a vector space definition (see 
Sect. |II]) and the above classification of integrable matrices into types, which allowed us to derive 



ype 1 matrices, map them to Gaudin magnets and obtain an 



an explicit parameterization for 

exact solution for their spectra js]. Further, we demonstrated that n-dependent energy levels of 
Type 1 matrices have at least one and at most (iV — 1)(A^ — 2)/2 crossings, while levels of integrable 
matrices of higher types do not have to cross, even though they frequently do. 

In what follows we first lay out the general strategy for resolving the commutation relations 
in Sect. [Ill In Sect. IIIII we briefly review the Type 1 construction of Refs. and 19. Next, we 



obtain all Type 2 matrices in Sect. IIVI In Sect. |V] we construct integrable matrices of arbitrary 
Type M, which includes M = 1 and 2 as particular cases. Exact eigenvalues and eigenstates of 
Type M are derived in Sect. IVII Next, we present in Sect. IVIII 'gauge' redundancies in the ansatz 
parameterization. In Sect. IVIIII numerical methods for generating random commuting matrices 
and determining their type, are presented and, thereby, the number of parameters necessary to 
specify a Type M commuting family is determined. In this section we will also discuss whether 
the construction of Sect. |V] generates all commuting families. Next, in Sect. IIX) we look at the 
frequent level crossings that occur in Type M matrices and relate their number to the properties 
of the discriminant of the characteristic polynomial of those matrices. We analyze the typology 



of various Hamiltonian blocks of the Id Hubbard model in Sect. |Xl We conclude the main text 
with the summary and discussion of open questions. Finally, various proofs, details, and tangential 
results are relegated to appendices |X] through O 

II. RESOLVING THE COMMUTATION RELATIONS 

We start with primary definitions and general observations. First, we define Type M opera- 
tors. As outlined in the introduction, the task of identifying commuting operators linear in a real 
parameter u in matrix representation reduces to finding a certain number, independent N x N 
Hermitian matrices (n) = T* + uV^ such that 

[H\u),W{u)]=0 for all u audi, j = 0,...,K -1. (1) 

Equating terms at all orders of u, we have 

[T ,V^] = [T^ , V'] , [T ,T^] = [V' ,V^] = 0. (2) 

As discussed above, we also require that there be no u-independent symmetry common to all H^{u), 
i.e. there is no 7^ al such that [il.,H^{u)] = for all u and i, where / is an x identity 
matrix. When such symmetry is present, H^[u) acquire block-diagonal structure and the problem 
reduces to that of smaller matrices. Note also that every commuting family contains the trivial 
member 

H^{u) = {a + hu)I, (3) 

where a and b are arbitrary real numbers. 

In light of the restriction on the ti-dependence, being independent of W{u) with i = . . . ,K —1 
simply means linear independence, i.e. Yli ^iH^W) with real Oj is zero if and only if all Oj = 0. Note 
that due to the absence of n-independent symmetries linear independence of H^{u) is equivalent 
to that of or T* separately. Further, K is defined as the maximum number of independent 
commuting matrices in a given family. In other words, any H[u) = T + uV that commutes with 
all H^{u) is linearly dependent on those same H^{u), i.e. 

K-l 

H{u) = aiH\u). (4) 

1=0 

We see that a commuting family is a i^T-dimensional vector space Vk where H^{u) serve as basic 
vectors. There is a certain freedom in choosing a basis in the vector space Vk, which will be 
exploited in Sects. |TTI] and HVl 
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It turns out that a typical matrix of the form T + uV has no nontrivial independent commuting 
partners hnear in ?x Q]; it commutes only with H^{u) = {a+bu)I and itself. Commutation relations 
Eq. ([1]) impose strong constraints on the matrix elements of each H^{u) such that the probability 
of finding randomly generated matrices T + uV with nontrivial commuting partners is zero. For 
example, in the 3x3 case matrix elements of T and V must meet a single algebraic constraint for 
such a partner to exist [2^ and the number of constraints increases with the size of the matrix. In 
our experience, finding an operator with even one nontrivial commuting partner is quite difficult. A 
brute force approach to generating all such operators numerically (see Sect. IVIIip already becomes 
computationally prohibitive for > 6. 

Therefore, the minimum (and typical) number of independent H^{u) in the family is K = 2. 
The maximum possible number can be shown to be -fC = A^, i.e. K ranges from 2 to A^. 
Following Ref. |6| we term the commuting families with the maximum number A^ of if *(n) - Type 
1, those with A^ — 1 commuting operators - Type 2, etc. up to Type A^ — 1. An arbitrary Type 
M family is one with K = N — M + 1 independent commuting N x N matrices = + uV^, 
which consist of the trivial independent element H^{u) = (a + bu)I and K — I nontrivial ones. 

Thus, our task is to identify nontrivial solutions of commutation relations ([2]). An important 
observation, which considerably simplifies this task, is that we have freedom to choose bases in 
two distinct linear spaces without loss of generality. First, we can choose a basis in the "target" 
space - the Hilbert space T-Ln on which Hermitian operators T* and act. Indeed, commutation 
relations ([2|) are invariant with respect to transformations {T^,V^) — t- {U^T^UjU^V^U), where U 
is any u-independent unitary operator. Second, we have freedom to choose a basis in the vector 
space Vk of commuting operators, i.e. to make linear transformations H^{u) — t- a^^H^ [u) with 
any real, u-independent and non-degenerate matrix (flj). 

It is convenient to choose the basis in the target space to be the common eigenbasis of the 
mutually commuting matrices V^. In this basis, are diagonal matrices whose nonzero entries 
we will denote d^, i.e. 

diagonal of = [d\,d2, ■ ■ ■ , d]^) . (5) 

The commutation relation [T*,!/-'] = [T-', y*] in this basis reads (d^ — dl^)T^^ = {dl, — dm)Tl^. 
This implies that an auxiliary antisymmetric matrix S whose matrix elements are 

C _ km km ff'\ 

>^km - ,i _ri - j j K'^) 

does not depend on i, i.e. is the same for all members of the family. Moreover, it follows from 



Eq. dlD that 

T = W'+[V\S] (7) 

where is a diagonal matrix (the diagonal of r*)[3]. Its nonzero elements we denote /|, 

diagonal of = (/{,/^,...,/^). (8) 

Now the commutation relations [y*, V^] = and [T*, V^] = [T^ , V^] are automatically satisfied and 
Eqs. ([2]) reduce to a single equation [T^T-'] = or, equivalently, to 

[[V\S],[V^,S]\ + [[V\S],W^] - [[V^,S],W'] =0. (9) 

We have, for fixed i and j, a single matrix "master" equation involving iN + N{N — l)/2 variables 
- 4N corresponding to V\V^,W\W^ and N{N - l)/2 corresponding to 5 - and N{N - l)/2 
constraints. In addition to this master equation, we will find that the quantity 

^ijki = SijSjkSkiSii + SikSkiSijSji + SiiSijSjkSki-i (10) 

is an important one for determining matrix type. We will see that for Type 1, I'ijki = 0, and for 
Type 2, i^ijki = GxiXjXkXi where G is some real parameter and Xj, j = {1, . . . , A^} are related to 
the matrix elements of . 

Next, let us discuss the choice of basis in the vector space Vk of K commuting operators. Note 
from Eq. ([7|) that a linear transformation H^{u) — )• a^jH^{u) translates into the same transfor- 
mation on diagonal matrices and and does not affect the matrix S, i.e. S is independent of 
the choice of basis in Vk and in this sense is a universal characteristic of the commuting family. 
Moreover, we will see that for Types 1 and 2 the master Eq. ([9]) reduces to a single equation on 
matrix elements of S. Once S is fully parameterized, and can be determined from Eq. ([9]). 
For Type M > 2 families, we do not attempt to resolve Eq. ^ generally, but use the Type 1 and 
2 resolutions to formulate a working ansatz that generates such families. 

As mentioned above, linear independence of W {u) is equivalent to that of . Given Eq. ([5]), this 
means that K = N — M + 1 vectors are linearly independent. Then, via a linear transformation 

— )• ^ a^jV-' we can go to a canonical basis in the operator space Vk such that K — 1 = N — M 
of diagonal matrix elements d\ are zero for every i. This is the basis we will use for Type 1 and 
Type 2 families in Sects. IIIII and ITVl respectively. 

All Hermitian Type 1 matrices were explicitly parameterized in Ref. 0. In constructing higher 
types here we will restrict our analysis to real symmetric as opposed to general Hermitian matrices 
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for simplicity. Specifically, we will choose the defining parameters of Type M commuting families 
to be real and such that the resulting matrices are real symmetric. We note however that all 
other properties of Type M operators we construct in this paper ~ i.e. commutation relations ([2]) 
and exact solution for the spectra of W\u) (see below) - persist even when the parameters are 
unrestricted complex numbers. 



III. REVIEW OF TYPE 1 MATRICES 

Here we review the main results for Type 1 [3] ~ families of N x N commuting operators, 
linear in a parameter u that have a maximum number (N) of linearly independent members, which 
were also termed "maximal operators" in Ref. 0. As discussed in the previous section, we can 
go from a general operator basis to a convenient canonical basis by taking linear combinations of 
H^{u) = T* + uV^. Since there are linearly independent diagonal matrices in the case of 
Type 1, we can choose the basis so that all d\. in Eq. ([5|) are zero except one, i.e. V^^, = = 5ik 
for all i. 

Consider Eq. ([9]) in the Type 1 canonical basis. It turns out that the left hand side of this 
equation is an explicitly mostly empty matrix. Evaluating its nonzero matrix elements, we obtain 
after some algebra 



f]-fl = -^. (11) 



where /c are distinct and are the diagonal matrix elements of T*, see Eq. ([8]). Consistency 
requires that a triangular sum of such relations itself vanish, i.e. (/j — fl) + (fl — fl) + (// — /j) = 0, 
from which follows 

i^ijki = 0, (12) 
where I'ijki is defined in Eq. (jlOp . This is the main equation defining Type 1 matrices The key 



insight in Ref. 



y is the understanding that Eq. ([12]) necessarily implies the parameterization 

S,k = (13) 

Sj - Ek 



21 



with unrestricted real parameters 7j and £j characteristic of the commuting family (see also Ref. 
for an elegant alternative derivation of Eq. (jl3p with the help of Pliicker relations). 

Once the matrix S is known, it is straightforward to calculate and matrices T* with the help 
of Eqs. ([9]) and (fTTj) . A general Type 1 matrix H{u) = '^i^=idiW{u) with arbitrary real di takes 



the form 



Imln 



U dr, 



dm d'l. 



ni ^ n, 



dm, dj 



(14) 



Eq. (jl4p provides a complete explicit parameterization of all Type 1 matrices. Note also that the 
parameterization (I14p is independent of the choice of basis in the operator space Vk- 

In Ref. Ifl we further show that Type 1 operators map onto a sector of Gaudin magnets in the 
presence of external magnetic field and as a consequence of this correspondence, we use the well 
known exact solution of the Gaudin model^, S, to furnish an exact solution for Type 1. The 
components of eigenvector Vm{u) of Type 1 matrix H{u) are 



Xr 



(15) 



corresponding to eigenvalue 



TV 



Em{u) = 



fc=l 



47^ 



(16) 



where the Aj, i = 1, . . . , are subject to a single algebraic equation 

N 



E 



7| 



Xr, 



U. 



(17) 



An intriguing consequence of this exact solution is that a Type 1 H{u) typically multiply violates 



the Wigner^von Neumann non-crossing rule 
least once 



22 



28l | and, moreover, it is necessary that it do so at 



IV. COMPLETE PARAMETERIZATION OF TYPE 2 MATRICES 

In this section we construct all Type 2 families, i.e. linear vector spaces formed by — 1 
mutually commuting, independent matrices H^{u) = T* + uV^. The derivation is similar to that 
of Type 1 in the previous section, but somewhat more involved. We first choose a convenient 
canonical basis in the operator space Vk as discussed in the end of Sect. |TI] and evaluate matrix 
elements of Eq. ([9]) in this basis. This results in an equation constraining the matrix Ski defining 
Type 2 matrices, which we then solve. The final result - parameterization of Type 2 operators - 
is given by Eq. ([33]). 
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There are — 1 linearly independent diagonal matrices for any Type 2 family. Therefore, 
by taking linear combinations of W{u), we can go to a basis in the operator space where only two 
of the in Eq. ([5]) are nonzero. It is convenient to choose and parameterize these nonzero matrix 
elements as = l/xj for some fixed j and = — l/xj for i = 1, 2, . . . , j — 1, j + 1, . . . , A^. We also 
impose the following constraint on parameters x^'- 

N 

^Xfc = 0, (18) 

k=l 

SO that the identity can be formed by linear combinations of the W (u) , thereby ensuring that it is 
an element of the commuting family. 
Eq. ([9]) in this basis yields 

fi fi _ SjkSji SikSii . 

for k, I 7^ {i, j}. Similar to the analogue Eq. (jlip . consistency requires that a triangular sum of such 
terms vanish, i.e. {fl - //) + (// - /*J + (/^ - /^) = 0, from which fohows that XiVjkim = XjVikim- 
Without loss of generality we can define a quantity Gijki such that Vijki = GijkiXiXjX^xi, which 
implies that Gjkim = Gikim- This relation is general and must be true for all distinct i, j,k,l,m, 
and therefore consistently permuting indices yields Gij^i = const = G, where G is some nonzero 
real parameter. We therefore obtain 

t^ijki = GxiXjXkXi, (20) 

We will show that Eq. (j20p can be solved generically. Let us simplify by dividing out the x's, i.e. 
let Sij = Sij^fxi^Jx~j^ from which it follows that 



^ijki — Sij Sjj^Sf^iSii + SikSkiSijSji + SiiSijSji^Sj^i — G. (21) 



We are looking for a parameterization of Sij that resolves Eq. (j2ip . ideally in a manner that 
treats both of its indices similarly. Our strategy for pursuing this parameterization is to find 
another equation nontrivially related to Eq. (|21|) and use that to eliminate some of the S's. More 
specifically, Eq. (|2ip involves four indices whereas the S"s are two index objects - implying that 
Eq. (|2ip overdetermines the S"s. We are looking for an equation involving fewer indices from which 
to derive a parameterization. 

The Sij, like the Sij, are antisymmetric and, as it turns out, this property is sufficient to 
identically satisfy the equation: 

^mjklSijSikSii — VmiklSijSjkSji + l^mijlSikSjkSkl — ^mijkSilS jiS]fi = 0. 
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Given Eq. ()2ip . we can remove an overall common factor G such that 



SijSikSii SijSjkSji + SikSjkSki SiiSjiSki — 0. (22) 

Eq. (|22|) is now a nontrivial equation constraining the Sij beyond simple antisymmetry. We can 
use Eq. (j2ip and Eq. (j22p to eliminate S^;. This is the resulting equation 



^ N 2 ^ X 2 ^ x2 ^^^-l 

SijSiij + [SijSjij - iSiiSjij —G 



where we have chosen to separate terms involving index k and those involving index I. Note that 
both sides of the equation are identical but for the different indices, true for all k,l ^ {i,j}, and 
it is for this reason that we can define the object Pi j clS Si two-index quantity though it is defined 
through quantities involving three indices. From Eq. ()23p we obtain 



Sfj + Sjk) ' - [SikSjk + Pijf = G- p% , (24) 
for all k 7^ {i,j}- The key result of our strategy is that Eq. ()24p implies that 



Sik + Sjk = ~ cosh Xk, 

Sij 



(25) 



SikSjk = JG - sinhxfc - Pij, 



for some parameters Xk- 

Again, what we have done here is use the identically satisfied Eq. (j22p to reduce a four-index 
equation (|2ip to a number of three-index equations involving the index-less parameter G and the 
two-index object pij. Actually, at this point, there is a simple algorithm to completely determine 
all 5's. We can fix indices i and j and choose values for G, pij, Sij and these new parameters Xk, 
^ 7^ Given these chosen values, the Sik and Sjk can be determined using Eq. ([2^ . Once 

we have values for Sik and Sjk, substitution into Eq. (j22p allows us to completely determine Ski, 
k,l / i.e. 

n o SikSii - SjkSji E; cosh Xfc + Efc cosh 

^kl — —^ij ~ ~ ~ ~ — —, —, , (,^0j 

' S.kSjk-SuSji smhxk-smhxi 
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where 



r 



G-p: 



2 



cosh _4^jQ_ p2 gjj^]^ _^ 4^.^ ^ 



1 

and 5m = il accounts for an inherent ambiguity in sign. 

The above algorithm expHcitly favors two indices over the others, i.e. the fixed i and j above, 
and we would like to treat all indices in a unified way. Let us make the following change of variables: 
sinhxm = {q^m + f) / (sEm + t) for q, r, s, t such that qt — rs= 1 [29t], and m / {i, j}. By substituting 



this into Eq. ()26p and removing an overall scale (again without loss of generality), we find that 

Ski = Jx/(Ai-efc)(A2-efc)(Ai-eO(A2-£;) ^A±Ii, (27) 

given 



= r{em) = 6m\ 1 + , + , , (28) 

V Ai — Em A2 — Em 

where the dependence of Ai, A2, -Pi, P2 on parameters q, r, s, t, G, pij, Sij can be determined with a 
bit of algebra. We are not particularly interested in this functional dependence, however, because 
Eq. (f28]l generally satisfies Eq. (f2T]) . i.e. if we choose parameters Ai, A2, -Pi, P2 and Em for 1 < m < 
N, all Ski given by Eq. (|27l) . k,l = {1 . . . N}, satisfy Eq. (|2T]) where in these preferred parameters 
G = -PiP2(Ai - A2)Vl6. 

To define the Type 2 Sij = Sij ^Jxl^Jxj fully, we can simply choose .A^ — 1 of the Xj parameters 
and determine the remaining one so that Eq. (llSp is satisfied. We can, however, write Xj in a more 
transparent form. Let us introduce new real parameters 7^ such that 

(29) 



{\\- Ej){\2- Ej) 

By partial fraction decomposition 



(30) 



(Ai — Ej) (A2 — Ej) Ai — A2 I Ai — Ej X2 — £j I 

we see that we can impose Eq. (jlSp by requiring that Ai and A2 be among the N solutions to 

Ex^ = ^' (31) 



where B is an arbitrary real number and Em are understood to be distinct. Given Eq. (I27p we find 
that Type 2 S^i take the simple form 

Ski = l^^{Tk + ri). (32) 
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Fully parameterizing Type 2 Sm is the lions' share of the Type 2 parameterization effort and 
with a bit of algebraic manipulation we determine the fl , for all k using Eqs. Q and (|19|) . which 
are linear in Linear combinations of the resulting W[u) yield a general Type 2 matrix 



\H [u)\^^ = -fmln — 7, , m / n, 



(33) 



where 



7V-2 



dm = 50 + E (34) 

j = l J 

and are arbitrary real numbers. 

Eq. (jSSp together with auxiliary Eqs. (j28p . pip and p4p provide a complete parameterization 
of all Type 2 commuting families. 2A^ + 3 independent real parameters Pi, P2, ■ji and £i specify 
the commuting family; — 1 real numbers Qi further specify a particular matrix within the family. 
The choice of all these parameters is unconstrained, except simple explicit restrictions on Pi and 
P2 to ensure that Eq. (j28p produces real F^, which we discuss in detail in the next section for 
general Type M. 

The choice of parameters is not unique. For example, a uniform shift — )• + c yields the 
same family. The entire group of "gauge" transformations on the parameters that leave Type 
2 commuting families invariant is explored in Sect. IVIII Interestingly, it turns out that in the 
case of Type 2 there is an additional gauge freedom not present for other types, which allows an 
alternative, particularly transparent parameterization, see Eq. (j60p . 

V. TYPE M - THE ANSATZ PARAMETERIZATION 

Explicitly constructing the most general Type M > 3 commuting families is difficult if at all 
possible (see the discussion Sect. IXII for more on this). However, it turns out that the Type 2 
formulas of the previous section can be generalized to obtain some Type M families with arbitrary 
M. The idea is to use the same Eqs. (p2]l and (p3]l to generate the matrices, but to extend Eq. (|28|) 
by including more poles in its radicand. Specifically, we take 



F^ = F(e„), Via) = 5ia) 
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" p 



•v, — a 
j=i J 



N-M 

dm = go+Yl ir^ ' (36) 

where 5 (a) = ±1, Pj are real numbers and Am, rn = {!> • • • are solutions of Eq. (j3ip . 

In Appendix|A]we prove that Eq. (i33]l with Tm and dm given by Eqs. ([35]) and (f36]) . respectively, 
indeed produce Type M matrices as defined in Sec. |TT1 i.e. families ofK = N — M + 1 linearly 
independent mutually commuting matrices with no n-independent symmetry. Eqs. (|33p and (j35p 
contain 2N + M + 1 parameters characterizing the commuting family - 2N 7j's and Ei's, M Pj's and 
parameter B. As discussed in detail in Sect. lVIll there are certain gauge transformations, such as a 
uniform scaling of ji and £i or a uniform shift of £{, that leave the commuting family invariant. This 
allows to fix three degrees of freedom meaning that the number of parameters needed to uniquely 
specify a Type M commuting family produced by our ansatz is 2N + M — 2. In Sects. IVTIIl and IXll 
we argue based on numerical evidence and other considerations that a general Type M > 3 family 
is uniquely specified by 2N + 2M — 5 parameters. This implies that our construction misses M — 3 
parameters, i.e. it does not yield all commuting families for types M > 3. 

A natural way to choose a basis for the ansatz Type M commuting family is to define the 
K -1 = N - M nontrivial H^u) such that Qj = 6ij in Eq. ^ for 1 < i < iV - M. For z = 0, we 
take H^{u) to be proportional to the identity matrix, as given in Eq. In other words, 

H\u) =T' + uV' is given by Eqs. 1^ and §5^ with dk 4 = T ^- • (37) 

Ai+M — £k 

for i = 1, . . . , N — M. In particular, V"^ is a diagonal matrix and 

diagonal of V = (4,4, . . . , 4) = f T ^ , T ^ , • • • , T ^ 1 • (38) 

A general member of the commuting family is 

N-M 

H{u) = g^H\u), (39) 

i=0 

consistent with Eqs. ()33p and (j36p . where gi are arbitrary real coefficients. 

Especially interesting in the structure of the ansatz is the form of Eq. (|35p . The presence of 
these Tj's in Eq. (j33p is what essentially distinguishes Type M, M > 1, from the Type 1 operators 
reviewed in Sec. IIIH and therefore what allows for the systematic decrease in size of the commuting 
family, i.e. these T's see us go from maximal operator Ski = lkli/{£k — £i) to 

Ski = \^^{Tk + Vi). (40) 
2£k- £i 
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for Type M. The successive inclusion of poles in the radicand of Eq. (j35p . i.e. increasing numbers 
of Pi 7^ 0, incrementally decreases the size of the corresponding commuting family. In particular, 
two such poles yields precisely the Type 2 matrices of Sect. |IV] (see in particular Eq. ([28|) ). i.e. 
for M = 2, the ansatz is precisely the complete parameterization of Type 2 given in the previous 
section. 

For M = 1, the ansatz offers an alternate, but equivalent parameterization of Type 1 matrices 
to that given in Eqs. (jl4p . To see this, consider the off-diagonal elements of a general Type 1 
matrix given by the ansatz in light of the equation T^ — = Pi{em — en)/(Ai — em)(Ai — 

(^)]mn = Imlnl^ ^TT — m / n (41) 



we can re-express Eq. (|4ip in the Eq. (|14p form, i.e. 



if we have 

2 wTi + X 



£ 



3 



7j 



Pi yTj + z 

7i 



(42) 



where w,x,y,z are arbitrary real numbers satisfying wz — xy = 1. Note that this constitutes a 
threefold redundancy in the ansatz parameterization of Type 1 matrices. In Sect. IVIII we will see 
that this redundancy is in addition to another threefold parametric redundancy that all ansatz 
matrices share. Demonstrating the equivalence of the diagonal elements of maximal matrices as 
defined by Eqs. (jl4p to those of ansatz Type 1 matrices, as given by Eqs. (j33p . proceeds similarly. 

The above is sufficient to show that Type 1 matrices, as specified by the above ansatz, are 
precisely the Type 1 matrices of Sec. IIIII Note, however, that there is a subtlety in the definition of 
dm- In particular, the ansatz seemingly stipulates (i.e. through Eq. (j36p ) a restricted set of allowed 
V, even for Type 1, whereas the Eqs. (|14p take the V to be diagonal matrices with arbitrary real 
diagonal elements dm- In fact, these are equivalent statements because for M = 1 there are Noj 



and N in Eq. (j36p related via an N x N matrix. This matrix is similar to the Cauchy matrix 
and can be similarly shown to be nondegenerate. Therefore, for an arbitrary choice of dm there 
always exists the corresponding choice of gj and vice versa. 

Note that Eqs. (I14p . (I33p and (I35p . AI = 1, are explicitly identical when Pi = and all Tm = 1- 
However, the ansatz allows for some F's to be —1 when Pi = 0. So for an ansatz Type 1 matrix 
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with Pi = - which we might call a 'Type 0' matrix - if Fj = Tj, Sij has the form Eq. ()13p . 
If, however, Fj = 1 and Fj = —1, Sij vanishes. It follows that if we take a Type 1 commuting 
family and tune |Pi| ever smaller, while keeping all other parameters fixed, we find that when 
Pi = we are left with pseudo-Type 1 matrices with block diagonal structure and, consequently, 
ti-independent symmetry. Indeed, these 'Type 0' commuting families will be the direct sum of two 
different Type 1 families. 



Ensuring reality of ansatz matrices 



Next, we discuss the conditions on ansatz parameters such that the resulting matrices are real 
symmetric. For unrestricted complex parameters 7^, £{, Pi, B and gi the above construction, 
Eqs. (j33p and (j35p . produces complex symmetric matrices for which all properties of Type M 
matrices save for reality, but including the exact spectra derived in the next section hold. 

Moreover, it can be shown that there exist complex parameters such that the corresponding 
matrices are real Type M, and where the complexity cannot be removed simply by multiplication 
by some common factor. A block of the Id Hubbard model in Appendix ICl provides an interesting 
example of such a situation. However, the most general constraint on complex parameters 7^, £i, Pi, 
B and gi to ensure real Type M matrices appears to be rather non-trivial and we will not attempt 
to derive it here. Instead, let us choose the parameters B, 'yj,£j, Pi and gi to be real and constrain 
them so that Eqs. (|33p and (|35p always yield real matrices. In this case the only problematic 
component for producing real matrices is the Fj - we still need the radicand in Eq. ()35p to be 
positive. 

First, we rewrite the expression under the square root in Eq. ()35p as follows 

j=i - nj=i (Aj - cr) 

i.e. we traded M parameters Pj for new M parameters (pj. Note that Pj's can be expressed in 
terms of (pj^s by computing the residues at the corresponding poles 

^ nfc^j(Afc-Aj) 

Let us order ej so that £j < £j+i- Recall that Xj are solutions of Eq. ([3T]l . It can be shown (see 
e.g. Ref. that they are all real and lie between consecutive £j. We have < Aj < ej+i for B > 
and Ei-i < Xi < Ei for B < 0, where i = 1, . . . , N. If we now take (pj, j = 1, . . . , M, to be also real 
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and such that 

Sj < (hj < Xj for B > 

J J (44) 

Xj < (t)j < £j for S < 

all terms j^^rf^ > 0, which is sufficient (though not necessary) to ensure the positivity of the 
radicand in Eq. [35] and reality of the matrix elements. 

VI. EXACT SOLUTION 

In this section, we present the exact solution of the Type M matrices parameterized by the 
above ansatz, see Sec. |Vl The exact solution of the eigenvalue problem for Type M matrices 



1 



see also Sec. IIIII A 



proceeds in a manner reminiscent of that of the Type 1 operators of Ref 
particularly striking feature of this exact solution is that, though these operators have explicitly 
more structure than Type 1 operators, the exact solution is nevertheless essentially just as simple, 
in that the entire eigensolution derives directly from the solution of a single algebraic equation. 
This is in stark contrast to the situation typical with Bethe's ansatz where there are a large number 
of coupled algebraic equations to be solved to derive a small portion of the spectrum. 

We want to determine the exact solution H^[u) ■ v"{u) = E^{u)v^{u), where W{u) is a basic 
operator specified in Eq. ()37p . v"{u) and E^{u) are its eigenvectors and corresponding eigenvalues, 
respectively. The ansatz eigensolution is as valid for Type 1, as it is for all ansatz Type M and 
so the components of eigenvector v"{u) must bear some resemblance to that of our existing Type 
1 exact solution. Let us modify Eq. (jlSp somewhat and assume that the j^^ component of this 
eigenvector is given by 

= [^{n)], = \^ (r(<7) + r,) . (45) 

Expanding in terms of our parameters and using Eqs. (j37|) and (|4Up . we obtain, after a bit of 
algebraic manipulation 



V 



3 



^i,r(a)^+r^_g(r(a)-i) 

^ (J £rn. 2 



Note, the second term is a j-dependent coefficient on the j^^ putative eigenvector component. This 
term must vanish in order for Eq. (j45p to constitute an eigenvector. Thus, given an eigenvector of 
matrix W[u) specified by Eq. (j45p the corresponding eigenvalue is 
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where a is such that the fohowing algebraic equation 

^ }.^k_ (r(a) + r„,) - I {T{a) -l) = u, (47) 

is satisfied. We should note, however, that it remains to be proven that Eq. (j47p generically yields 
N distinct solutions given both sign choices for T{a), see Eq. ()35p . Nevertheless, numerical analysis 
has yielded no counterexample. 



E 



u 



FIG. 1. Numerical energy levels of a 12 x 12 Type 4 operator = Y^i=o9i^^{''^) with parameters 

B,^j,ej,Pj, and randomly chosen over a uniform distribution of real numbers. The eigenvalues exhibit 
many level crossings, frequently violating the Wigner-von Neumann non-crossing theorem. 



The behavior of the spectrum of a typical Type M operator as a function of the parameter u is 
particularly interesting. As distinct from the spectra of random real symmetric operators linear in 
a parameter, random Type M matrices - whose parameters B, jjjEj, Aj, Pi are chosen from random 
uniform distributions - exhibit frequent violations of the Wigner-von-Neumann non-crossing rule, 
see Fig. [H which states that eigenvalues of the same symmetry do not cross as a function of a single 
coupling parameter. 

It is interesting to note that numerical observation of individual H^{u) suggests that they 
undergo concurrent crossings of — 2 eigenvalues at some value of u, see Fig. [2l This is a 
surprising phenomenon that one would rightly expect from Type 2 matrices [s^ but, perhaps, not 
general Type M matrices. 
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FIG. 2. Numerical energy levels of a 12 x 12 Type 4 operator W{u), i = 3, with parameters B, "fjySj, Xi, Pi 
identical to those in Fig. [1] Note that the eigenvalues exhibit — 2 = 10 concurrent level crossings at a 
particular value of u. 



We can explain this phenomenon by considering the H^{u) eigenvalues. At a given u, and for 
all a satisfying Eq. ()47p we find, through some algebraic manipulation 



m 



Aj — Em <7 — Em 2 Aj — fJ 



U 



m 



(48) 



Thus, at 5 = + r^) we find that, for a ^ Aj, the ci'th eigenvalue of W{u) is zero. 

Naturally, not all the eigenvalues can vanish, otherwise H^{u) is trivial. Indeed, setting a = \i 
and u = u m. Eq. (j47p . we find that it is solved for both sign choices for r(cr), i.e. Eq. ()47p has 
a doubly degenerate solution (T-{u) = o"+(n) = Aj, such that r((T+) = — r((T_). Numerically, we 
find confirmation of this phenomenon wherein N — 2 c's correspond to zero eigenvalues, while the 
other two eigenvalues at u are non-zero. 

VII. GAUGE REDUNDANCY IN THE ANSATZ PARAMETERS 

Given a set of real parameters B,'yj,ej,Xi,Pi defining a Type M commuting family one might 
ask whether these parameters are unique, i.e. does there exist a set B,^j,ej,Xi,Pi yielding the 
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same Type M commuting family? 

On inspection it is apparent that parameters are, in fact, not unique, e.g. a uniform shift of 
the £ (and A) parameters leaves Eqs. (j33p and (j35p unchanged. Similarly, one can find a variety of 
uniform scaling schemes under which all matrix elements are invariant. 

More interesting, however, is the following linear fractional parameter transformation: 

£7 = — , A,- = , 7,- = — r ? , P. = -7^7 (49) 

where w,x,y,z are arbitrary real parameters subject to constraint wz — xy = 1, <; = —z/y and 
r(<;) is defined in Eq. (|35p . It can be shown then that the above transformation yields the corollary 
transformation: 

r, = j^^. (50) 

By direct computation it can be shown that the Sij and, therefore, all off-diagonal matrix 
elements of the members of the ansatz commuting families are invariant under the above. Now 
consider the action of this transformation on diagonal matrix elements of basic commuting matrix 
H\u) = + uV^ given by Eqs. dSZ]) and dM])- For we have 



We see that the effect of the transformation on diagonal matrix is a uniform scaling by (1 — 
<;~^Aj)^ and shift by a trace, ^^"""^(1 — <;~^Ai)/, i.e. 

r = (1 - ?-Ui)V^ + ?-!(! - ?-U.)/. (52) 

Using partial fraction decomposition, see Eq. (|30p . and Eq. (j35p it can be shown that 



{i-c;-^Ky^H\u) = w 



u 



(1 - ?-iA,)2 ' ^ 2 



J (Ai-ej)(?-ej)^ 



i.e. the gauge transformation preserves the matrix elements of the Type M matrices, up to a 
uniform rescaling, a shift in the coupling parameter u and the addition of a multiple of the identity 
matrix. 

Type 2 

With respect to the ansatz parameterization, in addition to the invariance under the above 
transformation. Type 2 matrices have a parametric redundancy not shared by Types M > 2. 
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Numerical work undertaken to "reverse engineer" ansatz parameters directly from the matrix 
elements of M > 2 commuting matrices produced by the ansatz (see Appendix [B|) yield back the 
input parameters, modulo the transformation (|49p . However, Type 2 matrices exhibit a 1-dim 
gauge freedom beyond that of Eqs. (09]). 

To see where this extra redundancy comes from, consider again Eq. (|27p . Let us now allow £m 
to be a smooth function of a parameter t such that 

~ _ 1 f {em{t),en{t)) + r {en{t),em{t)) 



£m{t) - £n{t) 



(53) 



where F (a, /3) = Sm\/ (Ai — a) (A2 — a) {(pi — j3) (02 — /3) and </>i, </>2 are defined by the equation 

(</>! - a) (<A2 - a) _ ^ ^ Pi ^ P2 



(Ai — a) (A2 — a) Ai — a A2 — a 

for all a. By direct computation it can be shown that Smn is i-independent if 



(Ai - e^{t)) (A2 - £m{t)) {<Pl - £^{t)) {(t)2 - £m{t)) , (54) 



dt 

i.e. £m{t) is the elliptic function of t corresponding to Eq. (|54p with initial conditions £m{^) = ^m- 
Thus Smn is invariant under the transformation £m{^) — ^ £m(i)- 

A general elliptic function - i.e. one defined by Eq. ()54p wherein Ai, A2, ^i, 02 are arbitrary 
complex numbers - is related to a specific Jacobi elliptic function through a linear fractional 
transformation. Concretely, 

_ wksn{{tj+t)/T)+X 

yksn{{t,+t)/T) + z' ^''^^ 

where w,x,y,z are complex numbers such that wz — xy = 1, and 

W\l + X _ I W\2 + X _ I W(j)l + X _ W(j)2 + X _ , , 

— : — K , — — —k , ; — K ■ — K, (.50) 

yAi + z y\2 + z y(pi + z y<p2 + z 

T = zb/cy^ {—y'^k? + w"^) (— y^A;"^ + w'^), and k is any one of the four roots of the equation 

' k-k-^ \ ' ^ (Ai-0i)(A2-02) , s 

k + k-^j (Ai-(/)2)(A2-0l)' ^ ' 

Elliptic functions have a number of beautiful transformation properties. Interestingly, it can be 
shown, using the well known 'angle' addition formulae involving Jacobi elliptic functions j32l] that 



sn(t/r + k/4) = Tj/T {-z/y) , (58) 

where T{(j) is defined as in Eq. psp and k is one of the characteristic periods of these doubly- 
periodic functions. This is a surprising intermingling of £j and Ej, the broader meaning of which 
remains mysterious. 
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The most compelling such exploitation of the behavior of Jacobi elliptic functions comes in 
reconsidering Eq. ([55]) . By directly substituting Eq. ([5^ into Eq. ([55]) we find that 

^ ^ 1 cn ((t^ + t)/r) dn + t)/r) + cn + t)/r) dn {{tm + t)/r) 

2r sn((t^ + t)/T)-sn((t„+t)/r) " ^ ^ 

However Smn is independent of elliptic parameter t. We can make this t-invariance manifest by 
using the addition law properties of Jacobi elliptic functions 3^ , whereby it can be shown that 



cn dn 



Smn = ^ /, \ (60) 

ZT / Lm Or, 



sn 



2r 



where sn^ + cn^ = 1 and k^sn^ + dn^ = 1. 

This parametric redundancy is peculiar to Type 2 - the naive generalization to higher types, 
i.e. extending Eqs. ([53]) and (|54l) by including Xi, 1 < i < M and generating ej(t)'s that satisfy 
hyperelliptic differential equations does not yield an invariant analogue of Eq. (j53p . That such a 
simple, geometric redundancy is embedded in the derived construction of Type 2 matrices is both 
beautiful and compelling, however its meaning is unclear. 



VIII. RANDOM GENERATION OF COMMUTING MATRICES AND COMPLETENESS 
OF THE ANSATZ 

In this section we address the issue of completeness of the ansatz ~ whether all Type M > 2, 
commuting matrices are also ansatz matrices, i.e. for any pair of commuting matrices both linear 
in a parameter n, do ansatz parameters B,jj,ej,\i,Pi and gi exist that yield back these matrices? 
As yet, we have no analytic means to directly answer this question, but numerically we find that 
the answer is 'Fes' for M = 3 and Wo' for M > 3. Specifically, numerical analysis shows that 
2N + 2M — 5 continuous parameters are needed to uniquely specify a generic Type M commuting 
family, i.e. the ansatz of Sect. |V]is short by M — 3 parameters. 

Below we will detail this analysis, wherein we 

1. generate two random commuting matrices and determine the type of the commuting family 
they belong to, 

2. then process one of these matrices through an algorithm designed to check whether a matrix 
belongs to an ansatz Type M commuting family and extract ansatz parameters if it does. 
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Recall from Sect. [TI] that any two commuting symmetric matrices H^{u) = + uV^ and 
H'^{u)=T^ + uV^ can be represented as {u) = uV^ + [V^, S] + and (u) = uV'^ + [F^ , S] + 
VF^, where ,V'^ ,W'^ are diagonal and S is an antisymmetric matrix. [H^ (u) , H"^ (u)] = is 
then equivalent to (see Eq. ([9])) 

[[V\S],[V'',S]] + [[V\S],W''] - [[V^S],W^]=0. (61) 

It turns out that we can generate a pool of generic commuting matrices by brute forcing a solution 



to Eq. (j6ip . Using Mathematica's FindRoot\3a\ function - an algorithm designed to search a D- 
dimensional parameter space for solutions to a set of D equations - we input 4:N random real values 
for the nonzero elements of V^,V'^, W^, W"^ and ask FindRoot to find a solution S to Eq. (j6ip whose 
matrix elements are close to those of a randomly generated antisymmetric seed matrix. Typically 
FindRoot is able to quickly find such an antisymmetric matrix. This solution appears to be one of 
many in a large discrete set of compatible antisymmetric matrices, not all of them real. We know 
that there are many solutions because changing the random seed matrix, given the same inputs 

,V'^ ,W'^ , frequently sees FindRoot land on an entirely different antisymmetric matrix. 

We believe the solution set discrete because when Mathematica's NSolve - an algorithm designed 
to search for many numerical solutions to an arbitrary set of constraints - attempts to solve Eq. ()6ip 
given those same AN inputs it manages to generate a discrete set of S. However if we decrease the 
number of inputs, e.g. specify only AN — 1 inputs, NSolve very quickly indicates that the system 
is under-constrained by at least one equation and randomly generates a linear constraint so as to 
proceed toward calculating a solution set. 

In practice, we observe that FindRoot algorithm can find at least one solution to Eq. ()6ip . 
given AN inputs. Notice, though, that this method of randomly generating commuting matrices is 
independent of matrix type, i.e. nowhere in the procedure does one explicitly specify the size of the 
commuting family to which H^{u) and H'^{u) belong. The initial 4iV parameters defining matrices 
V^,V'^, W^, W"^ and the antisymmetric solution matrix S correspond to some Type M commuting 
family, and we can determine the value of M by counting the number of linearly independent 
matrices in the commuting family. 

In particular we take a putative matrix H'^{u) = uV^ + + [^''^S'], where V'' and W'' are 
diagonal matrices whose nonzero entries, and respectively, are undetermined real variables, 
and solve commutation relation [H^{u), H''{u)] = 0. There will be A^(A^ — l)/2 equations linear 
in these 2N unknown and Let p be the number of these equations found to be linearly 
independent. Then it follows that, because H^{u) is a Type M matrix and, by definition, we can 
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only freely choose K = N — M + 1 of the non-zero matrix elements of and the overall trace of 
W^, p = 2N — {K + 1) , i.e. the commutator yields p = N + M — 2 linearly independent constraints 
on and W^. Thus counting the number p of independent equations tells us the matrix type, 
M = p- N + 2. 

We observe that this procedure for generating generic commuting matrices yields a distribution 
of types weighted toward the lower end, see Fig. [3l Interestingly, out of many hundreds of random 
generations, not one 6x6 Type 3 matrices was generated, even though we know that they exist 
and can be readily made by the ansatz. 



1 




6x6 7x7 8x8 9x9 10x10 

FIG. 3. The distribution of various matrix types 1 < Af < — 2 among randomly generated integrable 
matrices of sizes = 6, 7, 8, 9 and 10. AN random real values (flatly distributed between —1 and 1) are given 
to the nonzero matrix elements of diagonal matrices V^,V^, W-^, and Eq. (|6T1) is solved numerically. At 
each iteration of random generation, one such S solution is chosen among a large set of such, and matrices 
H^iu) and H^(u) are constructed and their type determined. We observe that random generation favors 
lower types, and we note with some surprise that there are types that appear with extremely low (albeit 
finite) frequency, e.g. 8x8 Type 6, or not at all, 6x6 Type 3. 

We cannot as yet prove that given some random diagonal matrices ,V'^ ,W'^ a solution 
S must always exist for Eq. ([6T|) . but for the reasons stated above we find empirically that 



25 



parameters are always sufficient and apparently necessary to generate two distinct commuting 
matrices. Accordingly, let us assume 4:N the correct number of parameters necessary to specify 
two distinct, generic N x N real symmetric commuting matrices, linear in parameter u, regardless 
of type. If those two matrices are Type M, of those 4:N inputs 2{N — M + 2) are used to uniquely 
determine the two independent members of the commuting family, i.e. 2(A^ — M + 1) for both 
unique V^s and an extra 2 inputs to fix a trace on both W^s not constrained by the commutation 
relations. If, in addition, we remove an overall scale on S matrix elements, that leaves 2A^ + 2M — 5 
of our initial AN input parameters that remain to uniquely specify a generic Type M commuting 
family. 

Are all Type M matrices also ansatz? 

We already know that all Type 1 and 2 matrices conform to the ansatz - this is essentially 
what the ansatz was designed to do, see Sects. IIVI and IVl Now let us consider ansatz Type M > 3 
matrices. Counting ansatz parameters B,^j,£j, Pi and removing the parametric redundancies (see 
Sect. IVIip we are left with 2A^ + M — 2 parameters to uniquely specify an ansatz Type M > 3 
family. Comparing to the number of parameters needed to specify a generic Type M commuting 
family, we see that ansatz Type M families appear smaller by M — 3 continuous parameters. 

These parameter counting arguments strongly suggest that a generic Type M commuting matrix 
is not an ansatz matrix, i.e. the ansatz generates only a small subset of Type M for M > 3. The 
case of Type 3 appears marginal, i.e. on the one hand the preceding argument indicates that there 
are sufficient parameters for ansatz to cover all of Type 3, on the other hand the ansatz is a naive 
extension on a full parameterization of Type 2 and there is no obvious reason for the ansatz to 
extend fully to Type 3. 

To gain further insight into the completeness of the ansatz, we randomly generate generic 
commuting matrices of arbitrary type (see above) and subject them to an algorithm (described in 
detail in Appendix [B]) that reverse engineers the ansatz parameters from the matrix elements if 
the matrix is ansatz Type M. However, when we run generic commuting matrices through this 
algorithm, all observed matrices of Type M > 3 are seen not to be ansatz matrices. In particular, 
in these cases the algorithm is unable to find nontrivial Tj satisfying Eq. ()B3p . 

What is surprising, though, is that the ansatz carries no guarantee that it covers Type 3, and 
yet every observed, randomly generated Type 3 matrix has an ansatz parameterization. In view of 
this result, we conjecture that the ansatz is complete for Types 1 and 2 - which we have proved - 
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and for Type 3, while incomplete for general Type M > 3 by an unknown set of M — 3 parameters, 
which we note empties for Type 3. In Sect. IXII we provider further arguments supporting this 
conjecture. 



IX. LEVEL CROSSINGS 

One of the signatures of parameter dependent integrable models is the presence of level crossings 
in their spectra. A general u-dependent Hamiltonian, H(u), obeys the Wigner-von Neumann non- 
crossing rule which states that when one plots the u-dependence of those energy levels sharing the 
same quantum numbers, the levels may approach one another closely, but they never cross 
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28l ]. The spectra of Type M Hamiltonians, however, behave very differently - they carry frequent 
(though not necessary [6 1) violations of the non-crossing rule, see Figs. [H [21 and[H 

For Type 1, we were able to offer a topological rationale for the necessary existence of at least 
one such crossing^], but the analysis in Type M is more difficult. We do observe, however that 
the number of crossings varies with type as 

(A^-l)(iV-2) ~ , , 

# of crossings = -M + 1-2L, (62) 

where L is some non-negative integer bounded such that the above expression is non-negative, and 
where M is an integer such that 

M > M. 

In the overwhelming number of observed commuting matrices M = M. Specifically, out of several 
hundreds of matrices analyzed, only two non-ansatz matrices were observed to be such that M ^ M. 
M then is overwhelmingly identical to type, but there are apparent exceptions (see Sect. lXlj ). Note, 
we always find M = M for ansatz matrices. 

The crossings' story gets more interesting when one extends the notion of crossings to "complex 
crossings" , by which we mean complex values of u wherein at least two of the eigenvalues of an 
operator H{u) are equal. More concretely, complex crossings occur at values of u corresponding to 
the roots of the discriminant of the characteristic polynomial 'P{X,u) of H{u). The discriminant 
is defined as Yli^j (Ei — Ej)'^, where Ei are the roots of V{X,u), i.e. it is polynomial in the 
coefficients of V{X, u) and, therefore, polynomial in u. The order of this polynomial can be shown 
to be N{N — 1), and, therefore, there are exactly N{N — 1) complex values of u (roots) where the 
discriminant vanishes, independent of type or even whether H{u) belongs to a commuting family. 
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However, from numerical analysis we see that for a general H(u) these roots do not correspond 
to eigenvalues being equal one another; they correspond to a situation in which there ceases to be 
a full eigen-decomposition of the operator. For a completely general, non-integrable matrix, the 
roots of the discriminant correspond to a spectral decomposition of H{u) involving Jordan blocks, 



rather than eigenspaces, see Ref. 



35l . A Type M matrix, W{u), however, is rather special. Its 



discriminant ~ which we denote Discj;^i(„) - factorizes as 

Discj^.(„) =P{u)Q\uf, 

where P{u) and Q^{u) are polynomials. It turns out that Q^{u) is a polynomial of order {N — 
1){N — 2)/2 — M + 1 and its roots are often real, but not always. Its complex roots come in 
complex conjugate pairs, hence the 2L term in Eq. (j62p . All of its roots correspond to values of u 
where at least two eigenvalues are equal. Moreover, if we consider two Type M matrices, H^{u) 
and H^{u) belonging to the same commuting family, i.e. [W [u) , H\u)] = 0, their discriminants 
will have distinct factors Q^{u) and Q^{u). However, they will have the same factor P{u) and its 
2A^ + 2M — 4 complex roots correspond to values of u where a breakdown of the eigen-structure 
to Jordan blocks occurs. 

Interestingly we also observe that if we express H^{u) as a polynomial in H^{u) - which we 
can do for any pair of non-degenerate commuting matrices, see Ref. - the coefficients of that 
polynomial are rational functions in u the denominators of which are all Q^{u), i.e. 

N-l jyj , . 
m=0 ^ ^ ' 

where Rm{u) are polynomials in u of order (A^ — 1)(A^ — 2)/2 — M — m + 2. Note that the failure of 
polynomials of a degenerate matrix to span that same matrix' commutant is reflected in the fact 
that the above polynomial's coefficients blow up at the roots of (5*(n), i.e. at those values of u 
where H^{u) undergoes a level crossing. 



X. TAXONOMY OF TYPES IN THE ID HUBBARD MODEL 

In Sect. IVIIII we explored the distribution of general and ansatz matrix types among randomly 
generated commuting matrices. Here we study the prevalence of types in a well-known integrable 
system. As mentioned in the introduction (Sect. U]), sectors of such systems corresponding to a 
complete set of quantum numbers are Type M operators. Our aim is to determine the values of 
M for different sectors and to see if these operators conform to the ansatz of Sect. |Vl 
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FIG. 4. The M-dependent spectrum of a randomly generated 7x7 Type 4 operator H{u). Generally the 
discriminant of the characteristic polynomial of an x matrix T + uV is a polynomial in u of order 
N{N — 1). We observe that discriminants of Type M matrices factor as Disc//(„) ^ P{u)Q{u)^ . While 
the complex roots of P{u) correspond to values of u for which the spectral decomposition of H{u) involves 
Jordan blocks, the crossings occur at roots of the {N — 1){N — 2) — A/ + 1 = 12 order polynomial Q{u). This 
Type 4 matrix has 10 crossings at u = {-3.37, -0.90, -0.65, -0.55, -0.52, -0.49, -0.36, 0.31, 2.37, 6.10} and 
we numerically confirm that they are roots of Q{u). Its other 2 roots are complex conjugate, corresponding 
to 'complex crossings' occurring at u = {0.61 + 1.95i, 0.61 — 1.95i} 



As an example, we consider Id Hubbard model on 6 sites, 3 spin up and 3 spin down electrons. 
Periodic boundary conditions are assumed. The Hamiltonian reads 

^ = r E E (sV^+i ^ + 4+1 + E f^it - ^) i^ji (64) 

where cj^, and Cjs are fermionic creation and annihilation operators, respectively, s denotes the 
spin projection, hjg = cj^Cj^ is the number operator and T and U are real parameters. There is 
a hierarchy of conserved currents commuting with each other and with the Hubbard Hamiltonian. 
The first nontrivial current is linear in both T and C/, while higher currents are polynomials in U 
of order 3 and higher. Following Ref. |34| . we choose energy units so that U — 4T = 1, which is 



equivalent to the replacement 

U = u r=(n-l)/4, 
where u is a dimensionless real parameter. 
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Hamiltonian (j64p for 3 spin up and 3 spin down electrons can be represented by a square matrix 
of size (C|) = 400. It has numerous parameter independent symmetries such as total momentum, 
spin, particle-hole symmetry etc. Using these symmetries, one can bring the Hamiltonian to a block 



2C and 



34 for details as well as the 



diagonal form with block sizes ranging from 1 to 16; see Refs. 
quantum numbers corresponding to block sizes shown in Table HI With the help of this procedure 
and algorithms described in Sect. I Villi and Appendix[Bl we determine the smallest possible blocks, 
their type and whether they conform to our ansatz of Sect.|Vl The results are displayed in Table [H 
The total momentum takes integer values P = 0, . . . , 5 in our notation. Spectra for momenta P 
and 6 — -P are identical due to a spatial reflection symmetry. For simplicity, we show the results 
only for momenta P = 1, 5 and P = 2, 4. 



Momenta P 


= 1,5 






Momenta P 


= 2,4 




# of identical blocks 


N X N 


Type 


Is it ansatz? 


# of identical blocks 


N X N 


Type 


Is it ansatz? 


1 


8x8 


Type 3 


Yes 




2 


12 X 12 


Type? 


No 


2 


3x3 


Type 1 


Yes 




1 


14 X 14 


Type 11 


No 


2 


16 X 16 


Type 12 


No 




2 


4x4 


Type 1 


Yes 


1 


14 X 14 


Type 3 


Yes 




2 


2x2 






2 


3x3 


Type 1 


Yes 




1 


16 X 16 


Type 6 


No 



TABLE I. Blocks of Id Hubbard model are matrices of various types, some of which are described by the 
ansatz of Sect.|Vl Each block corresponds to a complete set of parameter-independent symmetry quantum 
numbers. The number of nontrivial commuting partners is TV — 7\jf — 1, where M is the type. For blocks in 
the tables this number ranges from 1 to 10, where 10 corresponds to the 14 x 14 Type 3 block. 



We make the following observations: 

1. Recall that N x N Type M matrix has A'^ — M — 1 nontrivial commuting partners linear in u. 
We note that blocks in Table [J have 1 to 10 such partners. This implies that the Id Hubbard 
model has at least 10 conserved currents linear in u, while only one such current has been pre- 
viously identified. All 10 currents can be written in terms of fermionic creation/annihilation 
operators cj^,Cjs with the help of projectors onto the corresponding sectors. The resulting 
expressions however might turn out to be nonlocal and rather cumbersome. 

2. The exact solution of Sect. IVII applies to all ansatz blocks in Table HI Therefore, in all these 
sectors Bethe's Ansatz solution for the Id Hubbard model^, [s^ - a large number of coupled 
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nonlinear equations - reduces to single equation (j47p . It is interesting to see explicitly how 
such a simplification becomes possible. 

In Appendix O we explicitly write down one of the ansatz Hubbard blocks - the 8x8 Type 3 block 
with momentum P = 1. 

XI. DISCUSSION 

In this paper, we have defined quantum integrable systems, linear in a parameter n, to be 
operators for which there exist a number of mutually commuting operators similarly linear in u. 
We introduced a classification of these quantum integrable systems according to type, where we 
defined a Type M operator as belonging to a commuting family with K = N — M + 1 linearly 
independent members. For Types 1 and 2, we were able to resolve the corresponding commutation 
relations, Eq. ([9]), yielding a complete parameterization. For higher types, we extended the Type 
2 parameterization to create an ansatz, see Eqs. (|33p - (|36p . that parameterizes a subset of all Type 
M commuting matrices. Moreover these ansatz matrices are exactly solvable through a single 
algebraic expression, Eq. ()17|) . 

In addition to this analytic approach to resolving commutation relations, we also brute forced 
random numerical solutions. We observed that the number of parameters necessary to specify one 
of these randomly generated commuting families is 2N + 2M — 5, whereas the ansatz involves 
2N + M - 2, i.e. 

M — 3 = of general Type M parameters) — of ansatz Type M parameters) 

When we attempt to find ansatz parameters that can reproduce a randomly generated commuting 
family's matrix elements, we find that this is possible for Types 1, 2, and 3 only; Types M > 3 
do not appear to be ansatz matrices, in general. Interestingly, when attempting to find ansatz 
parameters for the blocks of the Id Hubbard model Hamiltonian, we find that a number of them 
are indeed ansatz. Additionally, we looked at the unusual frequent violations of the non-crossing 
rule in the spectra of Type M matrices and found that the number of (possibly complex) crossings 
of a Type M matrix is precisely the same as the order of a polynomial in u whose square is a factor 
of the discriminant of that matrix' characteristic polynomial. The other factor of that discriminant 
is polynomial common to all members of the commuting family, and its roots correspond to complex 
values of u for which all H{u) cease to have a full eigen-decomposition, whereby these matrices 
have nontrivial Jordan canonical form. 
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The ansatz constitutes a significant step toward a comprehensive theory of parameter dependent 
finite dimensional integrable models. However, there remains much to be understood about general 
Type M operators. For example, it is believed that there is some to-be-determined equivalence 
between the existence of parameter dependent conservation laws, and the existence of an exact 
solution. This relationship is on full display with ansatz Type M matrices. However we could, for 
the sake of curiosity, consider a Type — 1 "commuting family" by following the ansatz exactly 
as written, but taking M = N — \. Note that this a very degenerate family because the number 
of independent elements is -ftT = 2, i.e. a single nontrivial operator and the identity. At first 
glance, such matrices do not seem to be different from most parameter dependent matrices in 
that their only commuting partner is the identity. Nevertheless, these Type — 1 matrices are 
exactly solved by same equations (j46p and ()47p . This suggests that conservation laws linear in the 
parameter are not necessary for an exact solution, though they may be sufficient. Interestingly, we 
find numerically that all ansatz Type — 1 matrices have commuting currents quadratic in u. 
The exact relationship between these two properties remains elusive. 

Most interesting is the question of how the ansatz fails to completely parameterize all Type M 
commuting matrices; i.e. why ansatz Type M commuting families involve M — 3 fewer parameters 
than generic Type M. Parameter counting alone tells us that the ansatz parameterizes all Type 
3 matrices, which we have confirmed numerically. Nevertheless, nothing in the derivation of the 
ansatz necessitates this Type 3 parametric completeness: this result is a complete (albeit pleasant) 
surprise. But why this surprise coverage of Type 3? Why does the ansatz completely parameterize 
Type 2 with a redundancy that can be made manifest through elliptic functions? 

Let us start with the last question first and let it be the foundation of the following speculation: 
Type M commuting matrices live on compact Riemann surfaces of genus g = M — 1. From this 
point of view. Type 2 matrices live on tori which, when equipped with a complex structure, become 
Riemann surfaces of genus 1. Working backward. Type 1 lives on the Riemann sphere. Similarly 
Type 3 lives on a Riemann surface of genus 2, a complex 2-manifold with two handles, Type 4 a 
genus 3 surface with three handles, etc. 

What do we get for this idle speculation? First, the redundancy of Type 2 ansatz parameters 
detailed in Sect. IVHI involves elliptic functions. One way to view such functions is as doubly 
periodic functions over the complex plane, but another way to understand them is as the basic 
functions that live on a fundamental parallelogram, where identification of the parallelogram edges 
implies that these functions live on a complex torus. This alone is a middling justification for such 
speculation in that it is limited to Type 2. 
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The Riemann surface conjecture becomes broadly instructive when we further conjecture that 
ansatz Type M matrices hve on the hyperelliptic subset of Riemann surfaces of genus g = M — 1. It 
is a well known result of the study of Riemann surfaces that general surfaces of genus g have complex 
structures parameterized by 3(? — 3 complex parameters, whereas hyperelliptic Riemann surfaces 
are much simpler and have complex structures parameterized by 2(7 — 1 complex parameters 



If we speculate that the ansatz consistently maps Riemann surfaces of genus g = M — 1 to Type 
M commuting families, but that the ansatz is presently restricted to that subset of surfaces that 
are hyperelliptic, parameter counting tells us that we should be short g — 2 complex parameters, 
i.e. that same M — 3 parameter gap between randomly generated versus ansatz derived Type M 
commuting families. Moreover, it is well known that most Riemann surfaces are not hyperelliptic, 
but that all surfaces of genus 0, 1, and 2 are — much the same way that the ansatz does not cover 
most types, but is complete for Types 1, 2 and (apparently) 3. 

In future work^| we will show that the eigenvalues of the ansatz Type M matrices belong 
to the topologically restricted vector space of meromorphic functions on hyperelliptic Riemann 
surfaces of genus M — 1. Moreover, the components of the eigenvectors will also be given by 
meromorphic functions on these surfaces and we will attempt to generalize the manner in which 
these functions are determined to arbitrary, non-hyperelliptic Riemann surfaces. There we will see 
that the well-known Riemann-Roch theorem [38j] justifies the correspondence between a Riemann 



surface's genus and the commuting family's type. In particular this theorem will show that the 
more general genus- type relationship is given hy g > M — 1, where the typical situation satisfies the 
lower bound. From this point of view the fact that sometimes M ^ M (see Sect. lIXp is reflected 
in the fact that g = M — 1. 
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Appendix A: Does the ansatz really parameterize Type M? 

Here we prove that the ansatz of Sect. IVl indeed yields Type M families of commuting operators 
as defined in Sect. |TT1 Specifically, we need to show that: 1) H^{u) defined by Eq. (f37|) admit no 
ti-independent symmetry, 2) [H^ (u) , (u)] = for all i,j < N — M and 3) the linear space Vk 
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formed by W{u) has dimension K = N — M -\- \, i.e. the only matrices that commute with all 
W{u) are linear combinations thereof. 

1. Absence of u-independent symmetry. Let us assume such a symmetry ft ^ al exists, 
i.e. [il, H'' {u)] = for all u and i, then 

[n,v'] = [n,T] = 0. 

However, for ansatz Type M matrices we can choose a basis (j38p on Vk such that = l/(Aj —Sk) 
where, by premise, all are distinct. Consequently an ansatz Type M commuting family has a 
basis in which all diagonal matrices are non-degenerate. It follows then that the only il. that 
commutes with all these is itself a diagonal matrix. By directly computing the commutator we 
see that for to commute with all 

7m7n(w,n " (^n) I 1 = 0, 

where = f^mm which is only generally satisfied if 

for all m, n = {1, . . . , A^}, i.e. if is a multiple of the identity. 

2. Mutual commutativity. Using the relation 

(rS. - r|) (r? - .) . - E E (,._.„k£-\)(a.-.) - 

and a bit of algebraic manipulation we find that 



CijkllmPsPt{>^s - >^tY{£k - - gm)(efc - gm) . 

where 

Iklii^i - Xj) 



a 



ijkl 



4(A, - ek){Xj - ek){Xi - ei){Xj - q) (r^ + 1) (F; + 1) 
Recall that the Aj satisfy Eq. (j3ip . Given this, it can be shown, using partial fraction decomposition, 
that 

(Ail gm)(Ai2 gm) • • • (Ajj, gm) 



34 



provided ii,i2, ■ ■ ■ ,ik are distinct. All terms summed over the index m in Eq. ()Aip contain a 
factor of the form (jA2p for k = 2,3 and 4, where Zi,i2, ■ ■ ■ ,ik are by construction distinct. Thus 
[i7*(n), //-^(li)] = 0, for all i and j ranging from to — 1. 

3. Dimension of the vector space formed by W{u). Next, we show that there are no 
members of this commuting family that are linearly independent of W{u). Toward this end, let us 
consider a general matrix H{u) = T + uV and specifically the diagonal matrix V and its nonzero 
matrix elements dm. First, we note that if there exists an H{u) in the commuting family that is 
linearly independent of the H^{u) = T* + uV^ , i = 0,...,-fC — 1, but whose V is linearly dependent 
on the y*, then there exists a linear combination of H^{u) and H{u) that is both u-independent 
and in the commuting family - a situation precluded by the absence of u-independent symmetry 
demonstrated above. 

More generally, in Ref. 



19|, Shastry shows that for arbitrary (not necessarily ansatz) commuting 
matrices - independent of type - there exists a necessary (but not sufficient) set of constraints on 
dm such that H{u) commute with all elements of the commuting family. The constraint is 



dj fJ-jiki + dk fJ.k;ij + di fii-jk + ^ di Vijki 

ij^j,k,l 



0, 



(A3) 



where I'ijki and fJ.j-kl are quantities characteristic of the commuting family as a whole and not 
specific to any particular one of its members. Uijki is defined as in Eq. (fT2]l and is fj,j-ki defined as 
follows: 



l^j;kl 



f}-fi 

4-4 



f]-n 

d]-d\ 



, 02 C.2 
-t- JjkJjl 



^ ^ SjmSm,kSklSlj-^ 

m^j,k 



d^ — d^ 



4 - 4 



^ ^ SjkSklSlm,Si 



d^ - d' 



mj 



dl 
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We have used the matrix elements of H^{u) to define fij-ki above, but it turns out that fij-ki is 
independent of index i. Let us now consider these general quantities as determined by matrix 
elements given by Eqs. (j33p and (j35p . By direct computation it can be shown that, given the 
Sec. |V] ansatz for Type M matrices 



M M 

EE 

t>s s=l 



16 



7| 



{Xs - ej){Xt 



n 

q={k,l} 



{Xs -eq){Xt - Eq) 



T^ijkl 



M M 

EE 

t>s s=l 



PsPtjXs - Xtf 

16 



n 



7^ 



q={i,j,k,l} 

Now let us compute Eq. (|A3p expressing dm as = 
always be done because the N x N Cauchy matrix Asm = (As 



{Xs - eq){Xt - Eg) 



N 



=1 



for some Qs ■ Note this can 



is invertible 



2M- With dm so 
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defined, we directly compute Eq. ()A3P to determine the necessary Qs such that [H'^ {u) , H {u)] = 0. 
We find that 



M M N 



9s 



It 



t>s s=l i=l 

which requires that 



PsPtiXs - \th]lhf 



0, 



Pt 



P 



Pt 



* At — Ej ^ Xt — £k * Xf; — El 

ETsPs -s—y l^sP-s ^sPs 



(A4) 



As — £j ^ Xs — £k ^ Xs — £l 
where summations over s and t run from 1 to M and = gg "^^i /{Xg — £i)^ ■ Note that terms 



containing Qs with s > M cancel from Eq. (jASp by virtue of identity (|A2p . 

For the determinant to vanish, the rows in the matrix in Eq. ()A4p must be linearly dependent. 
By the uniqueness of partial fractional decomposition it follows that = r, i.e. the determinant 
vanishes if (and only if) Vg does not depend on s. Thus we require that 

, -1 



9s 



f N 

E 



\ (As - ei) 



A2 



for 1 < s < M. 



(A5) 



It turns out that the degree of freedom associated with r amounts to variation in an overall trace 
as it can be shown that 



N 



N 



^-E^ E 



It 



A2 



rB-\ 



(A6) 



independent of the index m. Thus, it follows that the only V allowed in an ansatz Type M family 
are constrained to be of the form 

K 



V = Y, ^i^' + dl, 



(A7) 



i=l 



where = gi+M — 9i+M and gi+M is defined by the right hand side of Eq. (|A5p with s = i + M, 
whereby the ansatz necessarily parameterizes matrices of exactly Type M. 



Appendix B: Inverse problem: determining parameters given an ansatz matrix 

Here we detail an algorithm that, given an arbitrary N x N symmetric matrix H{u) = T + uV, 
determines if it conforms to the ansatz of Sect. |V] and returns the ansatz parameters when it does. 
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The algorithm is based on the central observation that the difference between a Type 1 Sij, see 
Eq. and that of an ansatz Type M matrix, see Eq. (|4Up . is the factor rjij = (Ej + Tj)/2. If a 

Type M commuting family has an antisymmetric matrix S for which no such factor can be found, 
we know it does not conform to an ansatz parameterization. If, however, such a factor exists and 
it can be determined, then we know that 

= ^ = (Bl) 

Tjij £i £j 

for some "jj and ej and Sjj so defined satisfies the equation 

^ij^jk^kl^H + ^ik^kl^lj'^ji + ^il^lj^jk^ki = 0. (B2) 

Moreover the reverse is true, i.e. it can be shown that Eq. ()B2p implies Eq. ()Bip for some 7j and 
Ej. Consequently, if Eq. ()B2p can be solved for Tj, we can essentially strip an ansatz Sij of its 
factor r]ij, and the ansatz parameters jj and Ej can be determined. Without loss of generality 
we can multiply Eq. ()B2p by an overall factor of lQriijrjiki]iir]jk'r]jir]ki and look for solutions to 
resulting equation 

(Ej + Efe) {Tj + Ti) SijSjkSkiSii + (Ej + E;) (E^ + E^) SikSkiSijSji+ 

(Ej + Tj) (Efc + Ti) SiiSijSjkSki = (B3) 

This is a massively overdetermined set of N{N — 1){N — 2){N — 3)/4! equations, quadratic in 
the N terms Tj. We note that the equations are actually linear and homogenous in the terms 
Bij = EjEj and that it is possible to attempt solving for the Tj by first finding the minimal set of 
linearly independent linear equations in Bij and then solving for Tj. In numerical practice, however, 
Mathematica's NSolve - programmed to find all sets Tj that satisfy the iV(iV-l)(iV-2)(iV-3)/4! 
equations of Eq. (|B3p - quickly finds nontrivial values for the Tj directly from the overdetermined 
quadratic equations, if they exist. We find that for matrices with a generic, non-ansatz S, solving 
Eq. (|B3p yields only trivial solutions, e.g. Tj = 0, for all j but one. Eor ansatz matrices where 
M > 3, however, there are two equivalent, inversely related non-trivial solutions to Eq. ()B3p . i.e. 

Ti = {aD,,(3D7^} 

where a, (3 are arbitrary complex numbers. Note that having real Sij does not generally guarantee 
that the corresponding Ej are themselves real. Note also that, despite having proved that Type 1 
and 2 matrices always have an ansatz parameterization, one can just as well attempt to determine 
their E's using solutions to Eq. (|B3p . However, in both cases the Tj cannot be determined uniquely. 
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In the Type 1 case, NSolve generates a warning indicating that there are not enough constraints 
for it to express a solution set and that it must generate three additional constraints in order to 
proceed; from Sect. |V] we know that these additional constraints amount to fixing the threefold 
parametric redundancy unique to Type 1. In the case of Type 2, NSolve generates a warning that 
it must generate a single constraint to proceed; this corresponds to a single redundancy in the 
parameterization of Type 2 matrices involving elliptic functions, which we detailed in Sect. IVTIl 

If nontrivial solutions Tj exist, the above procedure determines their values. It follows from 
Eq. ()Bip that the resulting Sjj satisfy a corollary equation 

(7,S,-fc)-i - (7iS,fc)-' + ilAj)-^ = 0, (B4) 

for some 7j to be determined. This we do by fixing i and j and choosing arbitrary values for 7j 
and 7j such that we can use equation 

7fc ^ = Sjfc ((7jSifc)~^ - (7iSjfc)~^) . (B5) 

to determine 7^ for all k 7^ From here, determining Sk is straightforward, i.e. 

ek = £i- ^p^, for all, k / i (B6) 

where we have yet another arbitrary degree of freedom in our choice of Ei . These three parametric 
redundancies, unearthed by this algorithm, constitute an interesting gauge freedom in all ansatz 
Type M commuting families, see Sect. IVIII 

If and when the algorithm finds values 7j, ej and Tj for a putative ansatz Type M matrix H{u), 
proceeding to determine ansatz parameters Xj is a matter of determining whether there exists an 
H{u) = uV + r in the commuting family wherein the nonzero elements of diagonal matrix V are 
of the form 

1 



dj 



X-Sj 



for A to be determined by the algorithm, see Sect.|Vl Recall from Sect. I Villi that in determining the 
size of a random Type M commuting family, finding all matrices that commute with H{u) reduces 
to + M — 2 linearly independent equations linear in its commuting partner's diagonal elements. 
Of these equations, we know that because there are only N — M + 1 independent members of the 
commuting family, exactly M — 1 of these equations can be found that involve the N dj alone. 
Consequently there will be M — 1 constraints on A of the form 



j=i - 
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where coefficients Cj are determined by the matrix elements of H{u). Solving these simultaneous 



equations reduces to finding all A that satisfy M — 1 polynomials of order — 1. Generally, such 
simultaneous equations polynomial have no solution. These ones derived from ansatz matrices, 
however, have an M-element solutions set {Ai, A2, • • • Am}- To be ansatz parameters, each one 
must correspond to the same parameter B satisfying Eq. (|3ip . 

Finding the rest of the M parameters Pj is a matter of solving N overdetermined linear equations 



see Eq. (|35|) . If the algorithm can determine these P^, H{u) is ansatz. Failure to determine any 
of these parameters uniquely (up to the aforementioned gauge redundancy), or any inconsistency 
with respect to the ansatz equations indicates that the matrix is not ansatz. For example, it is 
possible that the algorithm could find some Tj and parameters ej,7j given some matrix, and yet 
fail to have a consistent ansatz parameter B. In practice, however, all matrices tested through this 
algorithm failed to yield ansatz parameters at the F-stage in so much as the algorithm could not 
find nontrivial F^ consistent with Eq. (jB3p . That is, if there are non-ansatz Type M matrices that 
satisfy Eq. (j32p for some Fj not satisfying Eq. (j35p . they appear to be rather rare in a random 
ensemble of Type M matrices. 



Appendix C: An example of a sector of the Id Hubbard model described by the ansatz 

Here we explicitly write down the 8x8, momentum P = 1 block of the Id Hubbard model, 
which is an ansatz Type 3 matrix, see Sect. Kl This block is of the form T + uV , where V is 
diagonal matrix with the following entries: 




M 




(CI) 
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and 
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-^/6 
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3 / 



(C2) 



Using the procedure outlined in Appendix [Bl we find the corresponding ansatz parameters 

£ = {-1,0,5,-5/2,5,-5/2,-1,0}, 



7 



^ 1 1 - i, 3\/6 + i\/6, -1 - 7«, 4 + 8i, -7 + i, -8 + 4i, 1 + i, \/6 - 3\/6f | , 



2^6 



and 



2 5e| - (1 - 37i) + 4 - 3i 
^' ~ 5e| - (1 + 37z) ffj + 4 + 3i ■ 
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